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^SJ ' Abstract. Herbertsmithite, ZnCu3(OH)6Cl2, is the prototype candidate for a spin hquid 

behavior on a geometricaUy perfect kagome lattice. Its discovery and the absence of any 

' I ' ' evidence for spin-freezing down to the lowest probed temperature to-date enable one to explore 

O ' the properties of kagome-related physics in an unprecedented temperature range. We review its 

f, , properties and discuss some open issues. A tentative comparison to models is also performed. 
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I . 1. Introduction 

'rt ■ 
^ , Herbertsmithite ZnCu3(OH)6Cl2, a rare mineral recently extracted and named after G.F. 

Q \ Herbert Smith (1872-1953) who discovered the mother atacamite family in 1906 [1], has been 

, ^, ' synthesized for the first time in 2005, in M.I.T. [2]. It has triggered a major renewal in the 

\ search of a quantum spin liquid on the kagome Heisenberg antiferromagnetic (KHA) lattice 

C^ ' and opened the possibility to directly expose theories to experiments for the KHA ground state 

^i . which has eluded any definitive conclusion for the last twenty years. Altogether with triangular 

^r^ \ compounds k(ET)2Cu2(CN)3 [3] and more recently, (ET)Me3Sb[Pd(dmit)2]2 [5) it combines 

^D ' low dimensionality, a quantum character associated with spins 1/2 and the absence of magnetic 

^^ , freezing which was coined as the "end to the drought of spin liquids" (SJ E]. Herbertsmithite 

t^^ ' is the first experimental example of a kagome lattice where no order has been found at any 

^~ ■ temperature well below J through all experimental techniques [71 [S]- This illustrates the power 

^-H ! of the combination of an enhancement of quantum fluctuations for S = 1/2 spins with the 

L| ' frustration of antiferromagnetic interactions on the loosely connected kagome lattice to stabilize 

•r^ . novel ground states of magnetic matter. 

r> I Most of the models point to an exotic ground state and two broad classes of spin liquids have 

C^ ' been explored theoretically, in particular for the kagome antiferromagnet |10) . (i) topological 

spin liquids which feature a spin gap, a topological order and fractional excitations |TT]; (ii) 

algebraic spin liquids with gapless excitations in all spin sectors |12j . 

Here, we partially review some of the most relevant experimental properties of 

Herbertsmithite, and complement some aspects of our recent review on this compound [6J. 

We tentatively discuss which model could be consistent with these findings within the context 

of well identified deviations to the ideal Heisenberg kagome antiferromagnet. 

2. Structure and magnetic couplings 

Herbertsmithite ZnCu3(OH)6Cl2 is the x = 1 end compound of the Zn-paratacamite family 
Zn2^Cu4_a;(OH)6Cl2 [HE]. It can be viewed as a double variant of the parent clinoatacamite 





Figure 1. Simplified structure of ideal 
Herbertsmithite where only Cu (kagome top 
and bottom layers), Zn (intermediate triangular 
layer) and O which bridge between two Cu in 
the kagome planes, have been represented. 



Figure 2. Local structure emphasizing 
the coupling path within a kagome layer 
and between one Cu kagome site and one 
Cu substituted on the Zn site, named 
interplane site in the text. 



compound {x = 0), situated at the other end of the family. Starting from the latter, a Jahn- 
Teller distorted S = 1/2 pyrochlore, the symmetry first relaxes from monoclinic (P2i/n) to 
rhombohedral {Kim) around x = 0.33, leading to a perfect kagome lattice in the a-h plane 
with isotropic planar interactions; then, in the c- elongated x > 0.33 pyrochlore structure, 
the magnetic bridge along the c-axis between a-h kagome planes is progressively suppressed by 
replacing the apical Cu^"^ by a diamagnetic Zn^"*". 

In Herbertsmithite, the 119° Cu-OH-Cu bond-angle yields a moderate planar super-exchange 
interaction J ~ 180(10) K in comparison with other cuprates such as High T^s where the Cu- 
0-Cu bonding is linear (J ~ 1000 — 2000 K). Note that the former value is consistent with 
that observed in cubanes where a similar OH bridge links adjacent coppers [13]. Both neutron 
low-T data on the x = Q compound and a simple Curie- Weiss analysis of the variation of the 
high-T susceptibility with x, suggest that the interaction between apical and planar Cu is weakly 
ferromagnetic and estimated to be J ~ O.IJ. This is also in line with Goodenough-Kanamori 
rules which predict a weak ferromagnetic exchange for a Cu-O-Cu angle of 96.7° and with the 
presence of a weak ferromagnetic component in the frozen phase of Zn depleted {x < 0.66) 
compounds. 

3. Categorizing non-ideality and the related impact on physical properties 

The dominant term in the Hamiltonian is certainly the Heisenberg exchange interaction between 
neighboring spins within the kagome layer (see above). Yet, frustrated antiferromagnets are 
known to yield non-spin liquid ground states under small perturbations. In the case of 
Herbertsmithite where both disorder and Dzyaloshinkii-Moriya (DM) anisotropy are at play, 
the absence of any kind of order or freezing seems to be quite surprisingly robust. 



3.1. Dzyaloshinkii-Moriya interactions 

In the absence of any macroscopic high-T susceptibility results published on single crystals, 
ESR proves to be the most appropriate technique to investigate the local spin anisotropy. The 
inclusion in the Hamiltonian of a Dzyaloshinkii-Moriya interaction Dij.Si x Sj is necessary 
since there is no inversion symmetry between two adjacent Cu. D has two components, one 
perpendicular to the planes, Dz, and one planar Dp which act differently on the susceptibility 
and the specific heat. The effects are quite non linear and the large value of Dp invoked quite 
early after the discovery of Herbertsmithite to explain the rapid increase of the susceptibility 
at low T [14[ [T5] is far beyond the limit set by the ESR results. The small measured value of 



Dp r^ O.OIJ and sizeable D^ ~ 0.04 — 0.08 J |161I17] yield a fairly small correction, even negative, 
to the susceptibility as compared to the Curie-like upturn. Given this, the simple difference with 
series expansion at temperatures as high as J/2 argues for an extra term dominated by -if not 
only due to- out of plane defects [18\ . Definitely, the latter is responsible for the drastic low-T 
upturn of the susceptibility. 

Studies of the phase diagram of a quantum Heisenberg kagome lattice with DM interactions 
show that an original spin liquid phase can survive, at variance with the classical case where 
any minute amount of DM anisotropy will lead to long-range order y[9j. For spins S = 1/2 
a quantum critical point indeed appears for D/J ~ 0.1 which separates the liquid and the 
Neel states |20t \2T\ 122) . A peculiar variation of the susceptibility x ~ J'~°'^^ at the QCP is 
expected 
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Figure 3. ESR spectra taken at 300 
and 10 K. The broad linewidth is due 
to a large anisotropy. The dashed 
line shows the anisotropy of the g- 
shift tensor, hence the minute deviation 
from a pure Heisenberg character of 
the spins, quite common for Cu'^"^ ions. 
In the inset the Dzyaloshinskii vectors 
are represented. The Cu and O are 
respectively represented by small and 
large full circles (red and blue online). 
Adapted from ref. [16]. 



3.2. Out-of-plane defects 

The existence of out of plane magnetic defects due to the substitution of Cu^^ S = 1/2 onto the 
Zn^"*" non magnetic triangular layer has been the subject of controversies for a few years but this 
issue seems now settled. Based on electrostatic Jahn- Teller effect arguments, Cu^"*" was indeed 
expected to preferentially sit in the distorted octahedral environment of the kagome planes. 
For a perfect stoichiometry of three Cu to one Zn, which seems to be achieved according to 
Inductively Coupled Plasma Atomic Emission Spectroscopy (ICP-AES) characterization, this 
would lead to a a perfect Cu^"*" kagome plane and a non-magnetic Zn^"*" site in between the 
planes J2j. There is now a consensus that this argument cannot be used as an all and nothing 
argument. 

• Low-T-macroscopic measurements can indeed be described consistently within a scenario 
where the Zn site is occupied by 15-20% Cu^+. The upturn of the susceptibility [8], the 
partial saturation of the high field magnetization measured up to 14 T at 1.7 K [23], the 
linear field dependence of a Schottky anomaly in the specific heat [24] are strong indications 
of a Zeeman splitting of quasi-free S = 1/2 moments. 

• The low energy magnetic spectral weight in inelastic neutron scattering shifts linearly with 
field [2S]. 

• There were long-standing evidences that neutron diffraction spectra from various groups 
were best refined using a 15-20% occupation of the Zn site by Cu, which was also 
acknowledged recently by the MIT group. Yet the shallow minimum of the R- factor was 
casting some doubt on this determination. An additional very convincing argument came 



from the X-ray anomalous scattering data on a single crystal which allow to better quantify 
the level of Zn and Cu present at each site, namely 15(2)% of Cu was consistently found on 
the inter-kagome site [26] . 
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Figure 4. (Color on line) Two scalings for the low-T magnetization, M, measured up to 14 
Teslas. The 0.2 K data deviate from the scaling. Left: the contribution from the kagome planes 
which remains linear at high field has been subtracted. The remaining part is associated with 
quasi- free spins and the Curie- Weiss like scaling which works well down to ~ 1 K, suggests 
some magnetic interaction between the spins of the order of 1 K. Right: following the idea of a 
quantum critical behavior developed in [23,an equally good scaling is obtained. Yet, the fitted 
part does correspond to the now well established contribution from Cu defects on the Zn site 
which rules out the adaptability of a DM induced quantum critical scenario M T~^'^^ [22J for 
these data. 



Given this, one has to wonder how these quasi-free spins might couple with the kagome 
planes and how they might impact on the physics. As stated in the previous section, a small 
ferromagnetic interaction at least one order of magnitude smaller than the in-plane interaction 
should be anticipated. Surprisingly, susceptibility data and magnetization experiments at 
low-T do not show any sign of such a ferromagnetic interaction but rather point at a weak 
antiferromagnetic effective coupling of the order of 1 K between these "lone" spins [23]. Since 
a simple dipolar coupling would yield too weak an interaction of 10 mK for 15% S = 1/2 
moments, one is led to assume that this effective AF interaction might reveal some indirect 
coupling through the Cu's of the kagome planes. Using a reciprocity argument, one might infer 
that the lone spins have little influence on the kagome planes physics down to ~ 1 K. 

One should as well note in passing that it has been proposed at odds to interprete this 
low-T behavior of the susceptibility as an evidence for a quantum critical behavior under the 
assumption that it could be safely attributed to the kagome plane physics [2T]. Given the 
arguments of the previous section - the low-T macroscopic susceptibility is dominated by the 
inter-plane channel defect-, the latter condition certainly does not apply here. Fig. 4 shows 
the comparison between a Curie- Weiss analysis of the susceptibility and the scaling proposed in 
ref. [27J. They equally work well for T > 1 K but fail to interprete the 0.2 K data which signals 
some different physics of the defects in the low-T range. The muon relaxation rate also displays 
some peculiar change around IK. Such an analysis in terms of quantum criticality therefore can 
be safely replaced by the more sound scenario where a weak effective interaction couples these 
moments. 



3.3. In-plane vacancies : an open issue? 

While 15-20% of Cu occupation of the inter-plane site was found, depending on the technique 
used and likely on the samples, ICP-AES analyses have always reliably supported the idea 
of a perfect global stoichiometry of 1.00:3.00 Zn/Cu ratio. This leads to the straightforward 
conclusion that intersite mixing occurs and a 5-7% amount of Zn^"^ occupies the kagome planes, 
then creating spin vacancies. 

This seems to be quite disputed by the recent X-ray anomalous scattering data [26], where 
an upper limit of l(3)%Zn on the kagome layer was found. On the other hand, pushing to the 
extreme limit of 5% error on the 1.00:3.00 Zn/Cu ratio from ICP-AES reported on the same 
sample, would yield around 4% Zn on the kagome sites if 15% of the Zn sites are occupied with 
copper. We see that we are at the limit of the consistency between the two sets of data taken 
on the same sample, which likely require further inspections. 

Two other arguments based on ^'^O NMR spectra give evidence for the existence of in-plane 
defects I29|. 



• Typical spectra are displayed in Fig. [5j From a detailed analysis of the marked singularities 
of the powder lineshape which positions are governed by the magnetic shift and quadrupolar 
effects, two different sites can be identified. The most intense line at high T, named 
main(M), is associated with O sites linked to two Cu sites in the kagome planes. The 
other line (D) which is very sharp at low-T can be consistently interpreted as coming from 
O sites linked with one Cu and one non-magnetic Zn in the kagome planes. It is indeed half 
shifted at high-T as compared to the (M) line and corresponds to 20% spectral weight, ie a 
5% Zn occupation of the kagome planes. Its sharpness at low T is an additional indication 
that it corresponds to very specific sites in the kagome planes, the susceptibility of which 
is not at all related to the Curie-like contribution of the out of plane defects. 

• A broadening of the spectrum of the main line (M) is observed at low T. Would the kagome 
plane be free of any defect, one could still argue that: 

(i) The DM anisotropy could be responsible of this broadening on a powder spectrum. In 
the absence of any calculated field orientation dependence of the susceptibility due to DM 
anisotropic interaction, we crudely used the calculated susceptibility for an isolated dimer 
in the presence of DM interaction performed for SrCu2(B03)2 j28) . Since O probes two 
adjacent sites, the non staggered susceptibility should give birth to a powder shape which 
maximum relative width would be of the order of •^ gfiBHD^/J^- This represents ~ 10""^ 
contribution to the shift, much too small. 

(ii) The (7-factor anisotropy [IB] could be responsible for some broadening. Its value < 17% 
is not enough to explain such a huge low-T broadening and that should scale with the shift 
which is not the case. 

(iii) Out-of plane defects could be responsible for such a width. An asymmetric broadening 
~ 1/T would be expected through a direct coupling of these defects to the planes, which 
is not the case. Such a coupling is also ruled out by the features of the (D) line reviewed 
above. 

One can therefore safely conclude that the susceptibility in the kagome planes is not 
homogeneous and the ^^O linewidth is a clear signature of the existence of a staggered 
response around spinless defects. 

Based on this conclusion the relative 4:1 intensities of the (M) to (D) line argues in favor of 
5% defects in the kagome planes for the investigated sample. These are likely due to 5% spin 
vacancies but one cannot of course discard that the (D) line might result from local structural 
in-plane distortions due to Cu substitutions on the Zn site. Another honest still open possibility 
is that the samples investigated in NMR and those studied in X-ray are somewhat different. 




6.6 6.7 6.8 
H (Tesia) 



7.0 



Figure 5. (Color online) ^^O 
NMR spectra from 175 K to 
0.47 K from ref. [29]. The 
shift reference is indicated by 
the vertical line. (M) and (D) 
refer to the two oxygen sites 
in the kagome planes as de- 
scribed in the text. Large 
full green circles indicate the 
position of the center of the 
line that would be expected 
from XmacTo- At 175 K, a de- 
tailed analysis was performed 
and two contributions are ev- 
ident from the displayed fits 
(see text and section 4). The 
hatched area represents the 
contribution from (D) oxygens 
sitting nearby one Zn defect in 
the kagome plane. 



4. Evidences for a gapless ground state 

Given the dominant contribution of out of plane defects in macroscopic measurements, only 
site-selective probes can give some insight into the physics of the kagome planes. There are 
certainly two ideal probes of the kagome physics, Cu and O nuclei. Unfortunately Cu has never 
been detected in the whole T-range and the detection of a signal which was reported in [30] is 
certainly encouraging but which Cu, either in-plane or out-of-plane, is detected is still an open 
question. O therefore appears to be the most secure probe. From a simple comparison to other 
probes at high-T, where the macroscopic susceptibility is dominated by the kagome planes, O 
is found to be one order of magnitude at least better coupled than CI [30] and two orders of 
magnitude better than implanted muons in ^SR experiments [32] . 



4.1. Kagome susceptibility from ^'^ O NMR J2^ \31f 

Tracking the position of the main line gives the shift, hence one can extract the hereafter called 
"intrinsic" susceptibility from 300 K~ 1.7 J down to 0.45 K~ J/400. A broad maximum is 
observed between the high-T Curie- Weiss regime and the low-T regime where the susceptibility 
monotonously decreases between J/3 and J/20 to finally level-off at a finite value. The behavior 
in the intermediate T-range can certainly be associated with a strengthening of the short- 
range antiferromagnetic correlations. Fig. [6] illustrates the huge difference between the intrinsic 
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Figure 6. (Color online) Thermal variation of the susceptibility of the kagonie planes measured 
from the shift of the main ^^O NMR line (red open symbols) compared to the macroscopic 
susceptibility ( green solid line). The red dashed line is a guide to the eye for the low-T shift. 



susceptibility of the kagome planes and the macroscopic susceptibility which is dominated at 
low T by the inter-plane defects. 

Using the DM calculations in a dimer model cited above gives a strong indication that the 
DM interaction is not in an obvious manner the origin of the non zero susceptibility at T = 
and points to a gapless ground state in Herbertsmithite . . . unless for some not obvious reason 
J should be rescaled to a much weaker effective value, e.g. 0.1 J. This calls for further realistic 
inclusions of the DM interaction into the kagome Heisenberg hamiltonian. 



4-2. Dynamical susceptibility 

Both NMR Ti relaxation measurements [29\ [50] and inelastic neutron scattering (INS) [SJ |21] 
give also a strong support to a gapless scenario, at least in the field and temperature ranges 
which have been probed. NMR Ti measurements taken using various probes agree well from 
30 K down to 1.2 K. ^^O data extend down to 0.47 K under an applied field of 6.6 Tesla where 
a gapless picture is found to hold with a power law behavior T^ ~ 2^o.7i{5) Although less 
sensitive, INS data taken down to 35 mK without any applied field consistently reveal that 
x"{^) increases at low T with an uncommon divergent uj~^''^^^' law. 

One should as well note that a very different T-behavior is observed in /iSR experiments. A 
tiny increase of the relaxation rate is found below 1 K, to reach a plateau under zero applied 
field f7]. Since the dipolar coupling of the muon to nearby moments is an important channel for 
relaxation, these dynamical effects might be associated with interplane Cu. 



5. A tentative comparison to models 

Finding the ground state of the kagome Heisenberg antiferromagnet has represented a "Holy 
Grail" in the field of frustrated quantum magnetism. Two scenarios have been in balance since 
the early days, that of a spin liquid, gapped or ungapped, where quantum fiuctuations melt 
away any form of magnetic order and that of a honeycomb valence bond crystal combining 
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Figure 7. (Color online) Semi-log plot of the relaxation rates obtained for all probe nuclei, ^^O, 
^^Cu and ^^Cl. Inset: Linear plot of T^ with a power law fit for the low T ^^O data. Adapted 
from ref. 



dimer bonds within a 36-sites unit cell. The winner is of course that of lowest energy. As an 
example, the Valence Bond Crystal (VBC) model was found to be the best in 2008 through 
dimer series expansions [,33j whereas in late 2010 a Density Matrix Renormalization Group 
(DMRG) approach strongly supports a spin liquid ground state with a gap, with a possible 
value of J/20 |3l|. 

On the experimental side, Herbertsmithite does not favor any of these scenarios. Neither 
a gap is observed much below J/100, nor singlets are detected. The central difficulty is to 
understand better how much the ground state is perturbed by the presence of in-plane defects 
and DM anisotropy. We have selected a few proposed scenarios and discuss below whether the 
experimental data might cope with them or not in the realistic experimental situation. 

• Dirac gapless spin liquid: Using a fermionic approach in the context of fractional excitations 
allows to derive a specific heat C ~ T^, a susceptibility x ^ T, and the NMR relaxation 
rate 1/Ti ~ T'' [351 [36] ■ Although the low-T shift roughly scales linearly with T below J/6, 
the unstability of such a spin liquid against DM anisotropy rather argues against such a 
scenario. 

• Valence bond crystal/glass: Nothing in the data collected on Herbertsmithite supports 
the idea of a simple VBC state [33], which would combine dimer bonds as proposed for 
the honeycomb 36 sites unit-cell. Including DM interactions certainly relieves the singlet 
ground state constrain and could lead to some agreement with the experimental data but 
some effective interaction ~ O.IJ should be at play to explain the susceptibility, which 
remains unclear. Interestingly, in presence of disorder, a valence bond glass (VBG) could 
be stabilized where disorder might restore some susceptibility [32j. In this VBG framework, 
Raman scattering should give a very characteristic response. The defects are also predicted 
to yield free spins not simply localized singlets. Yet the recent DMRG calculations seem to 
disqualify such a valence bond based state, see below. 



Gapped spin liquid: This is the present winner in energy as found in a recent DMRG 
approach but a gap is predicted of the order of J/20 far too high to be consistent with the 
present experimental findings. 

What about the message from exact diagonalizations (ED) ? Besides the increase of the size 
of samples which can be studied, the computing progress enables now to calculate a large 
number of levels in any S sector of clusters up to A^ = 36 spins. This opens the possibility 
of a detailed analysis of the magnetism of the cluster. Numerical data of various sizes were 
recently scrutinized and lead to the conclusion that the singlet-triplet gap could vanish 
at the thermodynamic limit although as early as 1998, it was thought to be of the order 
of J/20 [38^ I39j . The analysis of the minimum energy in each spin sector leads to the 
conclusion that the KHA could be gapless with a susceptibility consistent with the NMR 
results, although 4 times larger [40j. 

DM interactions and disorder can also be included in ED calculations |4H I42j. As in most 
models, a singlet is found to localize in the vicinity of a spin vacancy - a relief of frustration 
on the corresponding triangle - and DM anisotropy again restores a small susceptibility. 
The longer distance response is staggered. Both features reproduce qualitatively the NMR 
spectra which respectively signal a specific low-T response in the vicinity of the spin vacancy 
and a broadening of the line which cannot be ascribed only to DM interactions without 
defects. 



6. Conclusion 
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Figure 8. Plot of the various energy scales at play in Herbertsmithite. 



In this paper, we aimed at clarifying some features associated with the defects. Out of plane 
defects seem not to affect much the physics of the kagome plane but they certainly obscure the 
macroscopic response. Including DM interaction in models, controlling the in-plane defects in 
order to extrapolate to a perfect kagome lattice are certainly the way to go for joining theory 
and experiment. At the present stage, the major result is that Herbertsmithite is a gapless 
spin liquid under the field and temperature conditions used up to now in experiments. Once 
the prehistorical age of Herbertsmithite will be over, spinless vacancies might play the role of a 
smoking gun for the physics of these compounds. Certainly much work is needed to reach that 
stage and the first synthesis of millimetric crystals are encouraging news for future research in the 
field . . . and to the most pessimistic reader who would be tempted to argue that Herbertsmithite 
is too disordered to shed light on the KHA ground state, EDs give a rather contrasting view 
that many experimental results can already be reproduced at least at a qualitative level! We 
certainly are at the stage where the complexity of the compound is well appreciated and one 



can now expect some considerable step forward in the materials improvement and the physical 
discussions on this compound. 
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